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Abstract. We have computed a spherically symmetric model for the interaction of matter ejected during the 
outburst of a classical nova with the stellar wind from its optical component. This model is used to describe the 
intense X-ray outburst (the peak 3-20 keV flux was ~ 2 Crab) of the binary system CI Camelopardalis in 1998. 
According to our model, the stellar wind from the optical component heated by a strong shock wave produced 
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when matter is ejected from the white dwarf as the result of a thermonuclear explosion on its surface is the 
^ i ' emission source in the standard X-ray band. Comparison of the calculated and observed time dependences of the 

mean radiation temperature and luminosity of the binary system during its outburst has yielded very important 
£NJ ■ characteristics of the explosion. We have been able to measure the velocity of the ejected matter immediately after 

the onset of the explosion for the first time: it follows from our model that the ejected matter had a velocity of 
, ~ 2700 km/s even on 0.1-0.5 day after the outburst onset and it flew with such a velocity for the first 1-1.5 day 

under an external force, possibly, the radiation pressure from the white dwarf. Subsequently, the matter probably 
became transparent and began to decelerate. The time dependence of the mean radiation temperature at late 
expansion phases has allowed us to estimate the mass of the ejected matter, ~ 10 -7 — 10 -6 M©. The mass loss 
rate in the stellar wind required to explain the observed peak luminosity of the binary system during its outburst 
has been estimated to be M ~ (1 — 2) x 10 Mo/yr. 

OO ' 

Key words, classical novae - X-ray emission - numerical calculations 

1. Introduction At this evolutionary phase of a classical nova, the 

C$ ' on-going hydrogen burning on the surface of the white 

1.1. X-ray Emission During Classical Nova Outbursts dwarf temporarily turns it to the so-called "supersoft" X- 

Explosivc thermonuclear burning of matter on the surface ^ sourcc ( see ' e.g.,MacDonald et al. 1985; Kahabka and 

of a white dwarf is generally believed to be responsible for van den Heuvel 1997 ) " thc cntire surface of thc white 

the phenomenon of classical novae (see, e.g., the review by dwarf Produces a flux that is approximately equal to thc 

Gallagher and Starrfield (1978)). During a thermonuclear Eddington one. In an optically thick regime and at typi- 

runaway, a large mass accreted on the white dwarf surface cal white dwarf sizes > 10 cm > the latter corresponds to a 

(up to ~ 1O- 4 M , Prialnik and Kovetz 2005) is ejected blackbody temperature of about 20-50 eV. 

and accelerated to high velocities, - 1000 - 4000 km/s Classical and recurrent (i.e., erupting several times, 

(see, e.g., Prialnik 1986; Prialnik and Kovetz 2005). e.g., Nova RS Oph) novae can also be the sources of stan- 

At thc initial evolutionary phases of a nova outburst, a dard (~ 2 - 10 keV) and hard ( > 20 keV) X-ray emission, 

high mass loss rate of the white dwarf envelope causes the The most popular model for the generation of such cmis- 

effective photosphere of the envelope to increase in size up s ^ on * s ^ ne shock model. 

to 10 12 - 10 14 cm, depending on the energy band (see, e.g., In the external shock model, a dense envelope that rose 

Kato and Hachisu 1994). Subsequently, the photospheric from the white dwarf surface as a result of the explosion 

radius decreases over several months, because the envelope and that was accelerated to velocities U ~ 1000 - 4000 

gradually becomes optically thin, until the entire envelope km/ s produces a strong shock wave in the stellar wind 

becomes completely transparent to emission. from the companion star. The shock-heated wind matter 

is the source of X-ray emission up to energies 5-80 keV 

* E-mail: kate@iki.rssi.ru (hv ~ m p U 2 /2). For example, the described mechanism is 
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believed to be responsible for the generation of hard X-ray 
emission from the recurrent Nova RS Oph (see, e.g., Bode 
and Kahn 1985; O'Brien et al.1992; Sokoloski et al.2006). 

Evidently, the external shock model is naturally linked 
to the existence of a dense medium near the binary sys- 
tem; otherwise, the emission measure of the shock-heated 
plasma (EM ~ n 2 dV ) would not be enough to pro- 
vide the necessary (observed) X-ray luminosity. However, 
in most cases, the companions in classical novae are 
unevolved main-sequence stars (Bath and Shaviv 1978) 
whose weak stellar wind is unable to produce a dense 
medium near the binary system needed for the external 
shock model. 

In the internal shock model, the ejected envelope has 
a velocity gradient as a result of which the faster inner 
layers catch up with the slower outer ones to produce 
a shock wave inside the expanding envelope (Lloyd et 
al.1992; O'Brien et al.1994; Mukai and Ishida 2001). The 
internal shock model was successfully applied to such clas- 
sical novae as Nova Herculis 1991 (Lloyd et al.1992) and 
Nova Velorum 1999 (Mukai and Ishida 2001). 

The internal and external shock models are not mutu- 
ally exclusive and can make comparable contributions to 
the total X-ray emission from the nova under certain con- 
ditions. The most important factor here is the acceleration 
law of the ejected envelope (here, we discuss only the ini- 
tial evolutionary phase of the nova) and the density of the 
stellar wind from the companion star. It is possible that 
not all of the envelope matter is instantaneously acceler- 
ated to the velocities measured at later envelope expansion 
phases from the parameters of optical and infrared emis- 
sion lines (see, e.g., Hynes et al. 2002; Das et al. 2006). 
Therefore, the inner layers of the envelope can have higher 
outflow velocities and, hence, an internal shock that affects 
significantly the observational manifestations of the nova 
in X-rays can be generated inside the ejected envelope ir- 
respective of the efficiency of the external shock model. 

The envelope acceleration process is important not 
only in explaining the X-ray emission from the nova, but 
it is also fundamentally important in understanding the 
physical processes in the burning envelope of the white 
dwarf. Thus, for example, the enrichment of the enve- 
lope with white dwarf matter (mainly carbon and oxy- 
gen) is very important for the formation of a rapidly ex- 
panding envelope in actual nova explosions. The energy 
release in the white dwarf envelope through thermonu- 
clear burning computed in ID numerical simulations with 
convectivc motions described by the mixing-length theory 
allows high envelope expansion velocities to be obtained 
(see, e.g., Starrfield et al.1985; Prialnik 1986; Yaron et al., 
2005). However, more detailed 2D and 3D burning calcu- 
lations do not allow such high velocities to be obtained 
without the assumptions about greatly enhanced carbon 
and oxygen abundances in the envelope (Kercek et al.1998; 
Glasner et al. 2005, 2007). Measuring the envelope expan- 
sion law would be the most important observational test 
of the white dwarf envelope burning model. 



The expansion law of the envelope ejected from the 
white dwarf is exceptionally difficult to test, because this 
process takes place within the first hours or even min- 
utes after the onset of explosive thermonuclear burning, 
i.e., long before the maximum optical light of the nova, 
i.e., in fact, before its actual detection in astronomical ob- 
servations. The observations of hard X-ray emission from 
several classical novae a few days after their maximum 
light in the optical energy band are not enough for this. 

Only recently has a unique opportunity arisen to test 
the nova envelope expansion models using the X-ray emis- 
sion from one of the most unusual transients in the X-ray 
sky, XTE J0421+560/CI Cam, as an example. 

1.2. XTE J0421+560/CI Cam 

The source XTE J0421+560 (later identified with the well- 
known variable star CI Cam; Wagner and Starrfield 1998) 
was discovered as an X-ray transient on March 31, 1998 
(Smith et al. 1998). The flux from the source rose to ^2 
Crab in several hours. Subsequent observations of the sys- 
tem revealed a variable flux in both optical and radio 
bands (see, e.g., Clark et al. 2000). 

The properties of XTE J0421+560 differed strikingly 
from those of ordinary X-ray transients (see, e.g., Tanaka 
and Shibazaki 1996). Various emission lines typical of the 
emission from an optically thin plasma with a temperature 
of 1-10 kcV (Orr et al. 1998; Ueda et al. 1998; Revnivtsev 
et al. 1999) but absolutely atypical of the emission from 
bright X-ray transients were detected in its spectrum. 

It was also found that, in contrast to ordinary X-ray 
transients (accreting neutron stars and black holes), the 
flux from XTE J0421+560 was highly stable at all Fourier 
frequencies - apart from a general decline in flux, only an 
upper limit on the variability of the source's X-ray flux on 
time scales of ~20-1000 s could be obtained (Revnivtsev 
et al. 1999). 

Subsequent studies of the source's behavior at late evo- 
lutionary phases (Ishida et al.2004), in quiescence (Boirin 
et al.2002), and investigation of the evolution of the ra- 
dio emission from the 1998 outburst (Mioduszewski and 
Rupen 2004) suggested that the intense 1998 outburst was 
the explosion of a classical nova in the binary system CI 
Cam and the compact object in XTE J0421+560 is most 
likely an accreting white dwarf. 

For a classical nova, the X-ray emission from XTE 
J0421+560 during the 1998 outburst is extraordinarily in- 
tense. For comparison, the peak X-ray flux from the widely 
known recurrent Nova RS Ophiuchi (Sokoloski et al.2006; 
Bode et al. 2006) was a factor of 10-30 lower than that 
from CI Cam during its 1998 outburst. A strong stellar 
wind from the optical component in the binary system, 
a B[e] giant star (Hynes et al.2002; Barsukova et al.2006; 
Robinson et al.2002), is responsible for such a high lu- 
minosity of the CI Cam/XTE J0421+560 outburst (we 
will show this below). It is the presence of such a dense 
medium around the expanding envelope of the nova that 
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Fig. 1. Left: light curve of the system in the 3-20 keV energy band during the 1998 outburst of CI Cam. The dashed 
straight line indicates the dependence L x ~ t 3 . Right: time dependence of the mean temperature of the emitting 
matter during the outburst. The mean temperature was measured from the ratio of the 3-5 and 5-20 keV fluxes. The 
first temperature measurement on 0.1-0.5 day of the outburst was made using ASM data (from the ratio of the 3-5 
and 5-12 keV fluxes); the remaining measurements were made using PCA/RXTE data. The dashed line indicates the 
dependence T ~ t^ 7 . 



produced the necessary emission measure (EM ~ n dV ) 
of the shock-heated plasma. 

The assumption that the explosion of a classical 
nova is responsible for the giant X-ray outburst of XTE 
J0421+560/CI Cam in 1998, along with the high X-ray 
intensity of the transient (the peak flux from the source 
reached ^2 Crab), which allowed a large set of high- 
quality observational data covering almost the entire evo- 
lution of the transient's X-ray outburst to be obtained, 
provides an excellent opportunity to test the most impor- 
tant ingredients of the external shock model for classi- 
cal nova explosions. A unique feature of the set of XTE 
J0421+560/CI Cam observations from the standpoint of 
studying the evolution of the classical nova envelope ex- 
pansion is the transient's light curve obtained by the All- 
Sky Monitor (ASM) of the RXTE observatory half a day 
before the X-ray flux peak. It allows the motion of the 
expanding envelope to be "seen" starting from several 
hundred white dwarf radii (~ 10 11-12 cm). The evolution 
of the X-ray outburst of XTE J0421+560, which, on the 
whole, lasted about 10 days, was well covered by RXTE 
observations (Revnivtsev et al. 1999). This also makes it 
possible to compare the predictions of the shock models 
during the decline in X-ray flux. 

This paper is the first in a series of papers that inves- 
tigate the behavior of the X-ray emission from CI Cam 
using hydrodynamic simulations of the expansion of the 



white dwarf envelope and its interaction with the stellar 
wind from the companion star. 

The goal of this paper is to demonstrate that the main 
features of the X-ray outburst of CI Cam can be explained 
in terms of the model of a shock wave passing through 
the stellar wind during a classical nova explosion. We in- 
vestigated in detail the following important points of the 
explosion: 

— the envelope motion at the initial expansion stages af- 
ter the explosion; 

— the deceleration of the envelope as it interacts with the 
stellar wind; 

— the contribution from radiative cooling to the change 
in the temperature of the stellar wind from the com- 
panion passed through the shock wave. 

1.3. Basic Characteristics of the Emission from CI 
Cam During the 1998 Outburst 

During the giant outburst of CI Cam in 1998, the X-ray 
emission from this binary system had several characteris- 
tic features that, as we will show below, can be explained 
in terms of the model of emission generation as the shock 
wave produced by the expanding envelope of the white 
dwarf passes through the stellar wind from the optical 
companion. 
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Fig. 2. Plasma cooling functions in two energy bands ob- 
tained from the APEC model. The functions were reduced 
to the total particle number density. 

— At the rise phase of the light curve, the time depen- 
dence of the source's 3-20 keV luminosity was L x ~ i~ 3 
(Fig.l, left). 

— The peak luminosity of the source was £ pea fc(3 — 
2QkeV) ~ 3 x 10 37 \d/2kpc) 2 crg/s, where d is the dis- 
tance to the source. 

— The mean (effective) radiation temperature of the 
source was constant and equal to kT ~ 10 keV for 
2 days from the outburst onset (Fig.l, right). 

— After the second day from the outburst onset, the ra- 
diation temperature decreased approximately as T ~ 

£-0.7-0.6 

2. The Method Of Numerical Calculations 

2.1. The Computational Scheme 

In this paper, we used numerical hydrodynamic calcula- 
tions to derive the time dependences of the effective tem- 
perature and X-ray luminosity of a classical nova. 

In our calculations, we used a ID spherically symmet- 
ric code in Lagrangian coordinates with a staggered mesh 
(the cell radius, velocity, and mass are determined at the 
cell boundaries, while the density, pressure, and internal 
energy are determined at the cell mass centers). A more 
detailed description of this computational scheme can be 
found in Janka et al. (1993). 

We disregarded the gravitational attraction from the 
white dwarf, since the shock wave in our calculations is 
generated at a distance, 5 x 10 11 cm, at which the es- 
cape velocity from the white dwarf, ~230 km/s, is much 
lower than the envelope velocity, ~3000 km/s. We also 
disregarded the influence of the optical star on the shock 
propagation. 

2.2. Radiative Cooling and Thermal Instability 

In our calculations, we took into account the ra- 
diative cooling of the stellar wind matter from the 
optical companion passed through the shock wave. 
The bremsstrahlung of an optically thin plasma was 



assumed to be the main mechanism of the radia- 
tive losses.We calculated the rate of plasma en- 
ergy losses through radiation using the APEC model 
( |http://h ea-www. harvard. cdu/APEC/REF) by assuming 
that ionization equilibrium was established instantly in 
the heated matter. The derived cooling functions are 
shown in Fig. 2. For convenience, these functions were re- 
duced to the total particle number density. The cosmic 
abundances were taken from Grevesse and Sauval (1998), 
the mean molecular weight of the plasma particles at such 
abundances is /i = 0.61, and n/n e = 1.93, where n is the 
total number of particles. To take into account the radia- 
tive energy losses, we used the 0.01-50 keV energy band in 
our calculations, which is sufficient for the temperatures 
expected in our model. The 3-20 keV energy band was 
used to construct the model light curves. 

A simple estimate of the stellar wind density at which 
radiative cooling begins to affect significantly the ther- 
mal balance of the post-shock matter can be obtained 
by comparing the radiative cooling time r ra( j through 
bremsstrahlung alone with the characteristic times in our 
problem, several days. For our estimate, we will take 1 
day. The gas in our calculations was assumed to be ideal 
with the adiabatic index 7 = 5/3: 
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where T is the temperature of the emitting plasma, n\ is 
the total density of the stellar wind matter passed through 
the shock wave, and A is the plasma emissivity (cooling 
function). For the observed temperature, kT ~ 10 keV, 
the plasma emissivity is A = 6.45 x 10~ 24 erg cm 3 s~ 1 and 
a radiative cooling time of 1 day corresponds to an unper- 
turbed stellar wind density uq = rii/4 ~ 10 10 cm~ 3 . 

In the case where the radiative cooling time of the 
shock-heated matter is shorter than the characteristic time 
of its expansion (several days), an isobaric thermal insta- 
bility develops (see, e.g., Field 1965). In ID calculations, 
this leads to a number of nonphysical effects. However, 
since the matter cooled below a temperature of ~1 keV 
ceases to radiate in the energy band of interest to us (3-20 
keV) , its specific temperature does not affect the results of 
our calculations as long as the thermal instability devel- 
ops at constant pressure. Therefore, when the cell cooling 
time was less than 0.1 day (which corresponds to temper- 
atures < 0.9 keV at typical densities in our problem), we 
switched off the cooling in this cell and, at the same time, 
excluded its emission from our calculations of the observed 
X-ray luminosity and the mean temperature of the X-ray 
emission. 

Let us now show that our approximation of the instant 
establishment of ionization equilibrium downstream of the 
shock is valid. Before estimating the time it takes for ion- 
ization equilibrium to be established, we must calculate 
the time it takes for a Maxwellian ion velocity distribution 
to be established and the equipartition time of the post- 
shock ion and electron temperatures r ie , since the shock 
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wave heats only the ions. Spitzer et al.(1965) showed that 
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where T$ and T e are the ion and electron temperatures 
in kelvins, is the ion density, A e = 1/1836 is the atomic 
weight of the electron, Ai is the atomic weight of the ion, 
and InA is the Coulomb logarithm. 

For our estimates, we will take Ai = 1, ^nA = 15, Tj = 
10fceV,T e = leV, = 1, and = n e . Then, 



95n i,iV 



Ti e ~ 0.04n ijl0 s 

where n^io = rii/10 w cm~ 3 . The time it takes for the 
electron temperature to be established is r ee = \[A~ e Tn = 
2-2^ s. 

To estimate the time it takes for ionization equilibrium 
to be established downstream of the shock, we will use a 
formula from Masai (1994): 
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where n e is the electron density. It follows from our 
estimates that the longest time it takes for the tempera- 
ture to be established (ja ~ 50s) and the time it takes 
for ionization equilibrium to be established (r eg ~ 50s) 
downstream of the shock are much shorter than the ra- 
diative cooling time of the matter (r ra d = 1 day) at the 



same density (here, rii 



m/2 = 2 x 10 1 
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Hence, our approximation of the instant establishment of 
ionization equilibrium downstream of the shock is valid. 
We also see from our estimates that the shock at matter 
parameters typical of our problem is not radiative (since 
Ta << T rac i), i.e., the density jump at the shock is defined 
by the formula n\ = ^jtiq and is ni/n = 4 for our value 
of 7 = 5/3. 



2.3. Calculating the Mean Temperature of the X-ray 
Emission 

Below, wc will show that, in general, the post-shock stellar 
wind matter is a multi-temperature plasma. Consequently, 
the radiation temperature that we measure based on X- 
ray observations is an average and it may not be equal 
to the temperature at the shock front. Tthcrcfore, to ob- 
tain the calculated mean temperature, we used the same 
averaging procedure as that for the observations. We cal- 
culated the ratio of the 3-5 and 5-20 keV fluxes, which, in 
turn, corresponds to a certain temperature in the radia- 
tion model of a single-temperature optically thin plasma. 
The flux ratio in the model of a singlc-tcmpcrature plasma 
was calculated using the APEC model. The efficiency of 
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Fig. 3. Spectrum of the matter passed through the shock 
wave composed of the spectra of the individual cells within 
which the matter has the same temperature (solid line). 
The dashed line indicates the spectrum of a single tem- 
perature plasma with the temperature determined from 
the flux ratio (see the text). 



the method is shown in Fig. 3, where the solid line in- 
dicates the combined spectrum of the post-shock matter 
at some instant of time and the dashed line indicates the 
spectrum of a single temperature plasma with the tem- 
perature determined by the method described above. We 
see from Fig. 3 that the method is efficient in our case. 

2.4. Initial Conditions 

The stellar wind density, temperature, and velocity pro- 
files were specified in our calculations as the initial condi- 
tions. The first inner cell of the computational grid with 
a finite mass was specified as the piston- envelope. This 
allowed us to subsequently consider the influence of the 
piston envelope mass on the shock dynamics. 

In our calculations with the piston-envelope, the ini- 
tial cell size was Ar = 5 x 10 11 cm. The inner bound- 
ary of the computational grid was placed at 10 10 cm. We 
performed calculations with various initial grid cell sizes 
to ascertain the optimal computational time/accuracy ra- 
tio. An increase in the resolution by a factor of 10 caused 
a change in the mean temperature of the emitting matter 
by a few percent while increasing significantly the compu- 
tational time. 

In our case, the stellar wind temperature does not af- 
fect the shock generation; therefore, it was taken to be 
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Fig. 4. Relative positions of the stars during the outburst. 
The hatched region corresponds to the possible positions 
of the white dwarf during the outburst. The black cir- 
cle marks the optical companion. The gray rings indicate 
the regions with a dense stellar wind; the wind density at 
the outer boundary of each ring is half that at the inner 
boundary. Arrow points the direction to the observer. 



approximately equal to the effective temperature of the 
optical star, - 10000 - 200001-iT : T = 1 eV. The temper- 
ature of the piston-envelope was taken to be equal to that 
of the stellar wind. 

Since the stellar wind velocity (35 km/s; Robinson et 
al.2002) is low compared to the envelope velocity (~1000- 
4000 km/s) and does not affect the shock generation ei- 
ther, it was set equal to zero in our calculations for sim- 
plicity. 

For the models with a constant envelope expansion ve- 
locity, constant velocities were maintained at both bound- 
aries of the envelope cell. 

When the shock wave was modeled in the Sedov 
regime (a strong point explosion in a spherically symmet- 
ric medium), the internal energy of the first four cells in- 
creased to the necessary explosion energy, whereupon the 
system relaxed. For our calculations in this regime, we set 
the initial cell size equal to Ar = 2x 10 11 cm. 

3. Positions Of The Components In CI Cam At 
The Time Of Its Outburst 

The optical component in the binary system CI Cam is an 
early B[e] star (Hynes et al.2002; Barsukova et al.2006), 
presumably a B4 III-V giant (Barsukova et al. 2006). The 
radius of the optical component is = 3.6 — 5.9i? = 
2.5 x 10 11 -4.1 x 10 n cm. 
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Fig. 5. Equivalent column for photoabsorption as a func- 
tion of binary inclination and mass loss rate by the optical 
companion. 

If the star is a supergiant (see, e.g., Hynes et al. 2002; 
Robinson et al. 2002), then its radius is considerably larger 
and can reach 3 x 10 12 cm. The orbital parameters of the 
binary were determined by Barsukova et al. (2006): its 
orbital period is P or b = 19.41 ± 0.02 days, the semi-major 
axis of its orbit is asinz ~ 70i?© = 4.8 x 10 12 cm (here, 
i is the inclination of the binary), and its eccentricity is 
e = 0.62 ± 0.07. The distance to the binary system is 
not yet known accurately; in different papers, it varies 
between 1 and 10 kpc (Hynes et al.2002; Robinson et al. 
2002; Ishida et al.2004). Barsukova et al. (2006) found the 
distance to the binary to be within the range 1.1-1.9 kpc. 
In this paper, we took 2 kpc. 

Using the time of periastron passage T = 52198.5 
MJD (Barsukova et al. 2006) and the outburst onset time 
50903.33 MJD, we calculated the orbital phase of the bi- 
nary, 4>orb ~ 0.273, and the separation between its com- 
ponents, rs'mi ~ 6.6 x 10 12 cm, at the outburst time. The 
measurement errors of the orbital period and eccentricity 
lead to inaccuracies in determining the orbital phase and, 
accordingly, the separation between the components. In 
Fig. 4, the region of possible positions of the white dwarf 
at the outburst time is hatched and corresponds to the 
following parameters: rsini ~ (5.4 — 7.5) x 10 12 cm and 
4>orb = 0.2 — 0.34. The position of the binary relative to the 
line of sight was found by Barsukova et al. (2006) based 
on the shape of the radial velocity curve. 

4. Stellar Wind From The Optical Component 

4.1. Estimating the Photoabsorption in the Stellar 
Wind from the Optical Component 

Robinson et al.(2002) estimated the mass loss rate in the 
stellar wind from the optical component in CI Cam to be 
M > lO~ 6 MQ/yr. Using the stellar mass loss rate in the 
wind, we can estimate the stellar wind density near the 
white dwarf and compare it with the photoabsorption in 
the X-ray spectrum of XTE J0421+560. 

Significant photoabsorption was detected in the X-ray 
emission from CI Cam/XTE J0421+560 in various observ- 
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ing periods (Revnivtsev et al. 1999; Boirin et al. 2002). 
The emission from XTE J0421+560 in quiescence (Boirin 
et al. 2002) most likely results from mass accretion onto 
the white dwarf. The equivalent column for photoabsorp- 
tion in this spectrum may be the sum of the photoab- 
sorptions in the stellar wind within the binary system 
and in the accretion column in the immediate vicinity of 
the white dwarf. In this case, the large column density, 
N H L - 5 x 10 23 cm- 2 , obtained by Boirin et al. (2002) 
does not allow us to unequivocally attribute this to the 
contribution from the absorption precisely in the stellar 
wind. 

To estimate the stellar wind density, we can use the 
column for photoabsorption at the beginning or at the 
peak of the light curve for the 1998 outburst (NhL ~ 
5 x 10 22 cto -2 ). In this case, the X-ray emission emerges 
far from the white dwarf surface and, hence, it is in any 
case free from this additional photoabsorption. 

To calculate the dependence of the column density on 
the mass loss rate by the optical component, we assumed 
that the stellar wind from the optical star was spherically 
symmetric about the star, had a velocity v W i ne i = 35 km/s, 
and its density n depended on the distance to the stellar 
center R as 



M 



= 1.4 x 10 11 rr 2 2 « 3 - 5 1 M_ 7 , 



where fj, is the mean molecular weight of a single par- 
ticle, m p is the proton mass, t\2 — r/10 12 cm, M-i = 
M I 10 -7 M Q / yr , and U35 = u/35 km/s. The relative posi- 
tions of the star and the line of sight with respect to the 
binary system were determined above. 

In Fig. 5, the contours indicate several photoabsorp- 
tion columns for various orbital inclinations and mass loss 
rates by the optical component. We see from the figure 
that the outflow rate of the stellar wind from the optical 
component at the observed NjjL = 5x 10 22 cm -2 can lie in 
the range 4 x 10~ 8 < M < 2 x 10~ 6 M Q /yr, while N H L = 
5 x 10 23 cm -2 gives a constraint M > X 10~ 7 M Q /yr . These 
estimates are consistent with the values from Robinson et 
al. (2002). 
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Fig. 6. Averaged stellar wind density profiles: for sini = 1 
(solid line) and sini = 0.6 (dash-dotted line). The dashed 
curve indicates the law of decrease n oc r~ 2 . In this calcu- 
lation, the mass loss rate by the optical star in the wind 
is M = 4 X 10" 7 M Q /yr. 



To make the first simplest estimates of the density dis- 
tribution of the interstellar medium as a function of the 
distance from the white dwarf, we used assumptions about 
the stellar wind motion similar to those made in the previ- 
ous section and calculated the mean stellar wind density in 
the volumes between concentric spheres centered at the lo- 
cation of the white dwarf during the outburst. The radius 
of the optical component in our calculations was taken to 
be 5.9-Rq. The derived average stellar wind density profile 
as a function of the distance from the white dwarf is shown 
in Fig. 6. We see that the density is almost constant near 
the white dwarf and behaves as a typical stellar wind den- 
sity profile with n oc r~ 2 starting from some distance. We 
will use this approximation of the density profile in our 
calculations. It also follows from the derived profile that 
r c shu ~ 6.6 x 10 12 cm. Below, r c denotes the distance at 
which a constant stellar wind density profile transforms 
into a decreasing one. 

In Fig. 7, the contours show how the stellar wind den- 
sity near the white dwarf (at r < r c ) depends on the 
orbital inclination and the mass loss rate by the optical 
component. 



4.2. Estimating the Density Distribution of the Stellar 
Wind Relative to the White Dwarf 

The stellar wind density at large distances r from the 
star generally decreases as n oc r~ 2 . However, since at 
the outburst onset the envelope that rose from the white 
dwarf generates a shock wave around the latter even inside 
the binary system, we cannot assume the density profile 
n oc r~ 2 (where r is the distance from the white dwarf) 
to be correct from the very onset of the outburst. 

At small distances from the white dwarf, r < r c , the 
stellar wind density will not decreases as r~ 2 (where r 
is the distance from the white dwarf) but will be more 
uniform. At larger distances, r > r c , the transition to a 
decreasing density n oc r~ 2 will occur. 



5. The Scheme Of Shock Waves During The 
Outburst 

The external shock model that we consider here suggests 
that the white dwarf envelope with a mass M e j moves 
under the action of an external force (possibly, the radi- 
ation pressure from the white dwarf on which thermonu- 
clear burning of the remaining accreted matter continues) 
for some time At after the onset of its expansion; subse- 
quently, it flics ballistically, i.e., the law of motion of the 
envelope is determined only by its interaction with the sur- 
rounding matter and is not the result of its acceleration 
by the pressure from the burning part of the envelope. 

As a result of the envelope expansion, an external 
shock travels through the stellar wind, while either a shock 
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Fig. 8. Scheme of shock waves generated by the expanding 
envelope: U is the velocity of the contact discontinuity or 
the envelope and D is the velocity of the forward shock. 
The subscripts "0" and "1" denote the quantities in the 
unperturbed stellar wind and downstream of the forward 
shock, respectively. 



travels through the envelope inward or a rarefaction wave 
is initially generated which subsequently transforms into 
a shock as the envelope expands and as the pressure in it 
decreases, depending on the envelope velocity for a given 
ratio of the initial pressures. 

A schematic view of the system of shocks is shown 
in Fig. 8. In principle, the X-ray emission can come not 
only from the forward shock but also from the reverse 
one. However, since the reverse shock should be weaker 
than the forward one and since the radiative cooling of 
the matter downstream of the reverse shock is strong due 
to its high density, the expected temperature of the matter 
downstream of the reverse shock will not be high enough 



for it to radiate in the X-ray band under consideration (3- 
20 keV, corresponding to temperatures T ~ 10 7 — 10 s K). 
Therefore, we assume here that the contribution from the 
envelope to the observed X-ray emission is negligible and 
disregard it. In our calculations, the envelope was consid- 
ered as a piston of finite mass producing a shock in the 
stellar wind and the letter U denotes the velocity of the 
piston-envelope. 

6. The Effect Of Nonuniformities In The Stellar 
Wind Density Profile On The Mean 
Temperature And Luminosity Of The Emitting 
Matter 

The properties of the X-ray emission from CI Cam dur- 
ing its outburst strongly suggest that it is formed by the 
emission of an ordinary hot optically thin plasma (Ueda 
et al.1998; Orr et al.1998; Revnivtsev et al.1999). The lu- 
minosity of a unit volume of optically thin hot plasma 
is proportional to the square of its density and, hence, 
is very sensitive to density variations around the white 
dwarf. In the binary system CI Cam, which consists of 
a young star with a strong stellar wind and an accreting 
white dwarf, there must exist various deviations from a 
spherically symmetric distribution of stellar wind matter 
around the white dwarf, such as the density jumps in the 
wake produced when the white dwarf moves through the 
stellar wind (Dumm et al.2000), the nonspherically sym- 
metric stellar wind from the optical star, etc. This was 
clearly demonstrated in a recent paper devoted to numer- 
ical simulations of the matter distribution in the recurrent 
Nova RS Oph, which also contains an optical star with a 
strong stellar wind (Walder et al. 2008). 

To estimate the effect of stellar wind density nonunifor- 
mities on the mean radiation temperature and luminosity, 
we calculated the envelope expansion with a constant ve- 
locity in a medium with a variable density that depended 
on the radius as n = 10° + 8 x 10 8 sin(ri 3 ), where r i3 
is the radius in 10 13 cm (see Fig. 9 ). For simplicity, the 
radiative losses of the post-shock matter were disregarded 
in the calculations. 
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The envelope motion with a constant velocity through 
a medium with a constant density produces a shock wave 
with a constant temperature at the front; consequently, 
the mean temperature of the postshock matter also re- 
mains constant (if the radiative losses are small). The lu- 
minosity of the shock-heated matter depends on the time 
as L oc t 3 (i.e., it is proportional to the volume of the emit- 
ting matter; a more detailed derivation of the formula will 
be given in the Section "Rise Phase of the Light Curve"). 

For a variable (changing by an order of magnitude) 
pre-shock matter density, the mean post-shock tempera- 
ture and the luminosity of the heated matter will behave 
in a different way (Figs. 10 (top left) and 10 (top right)). 
Figures 10 (bottom left) and 10 (bottom right) show de- 
viations of the mean temperature from the expected con- 
stant value and deviations of the luminosity from the ex- 
pected law L (x t 3 , respectively. We clearly see from Fig. 
10 that the mean radiation temperature is much less af- 
fected by density jumps: for deviations of the luminosity 
from the law L cx t 3 by a factor of 12, the radiation tem- 
perature changes only by a factor of 1.5. Therefore, to 
qualitatively describe the evolution of the X-ray emission 
from the 1998 outburst of CI Cam, we will primarily rely 
on the observed behavior of the emitting-plasma temper- 
ature. 

7. Mean Radiation Temperatures And Shock 
Velocity In The Initial Expansion Period 

Starting from the hrst hours of outburst development, the 
ASM/RXTE flux measurements for CI Cam in several 
channels allow the temperature of the emitting plasma 
to be estimated. We used the ASM/RXTE hard X-ray 
measurements of CI Cam in the 3-5 and 5-12 keV en- 
ergy channels, where the effect of photoabsorption clearly 
seen in the source's spectrum on 0.7 day (Revnivtsev et 
al. 1999) is weaker. Using the known hardness (flux ra- 
tio) of the Crab Nebula in the ASM 3-5 and 5-12 keV 
energy channels, we estimated what hardness of CI Cam 
corresponded to the radiation temperature of an optically 
thin plasma (the APEC model) within the first 0.5 day 
after the outburst onset. When modeling the radiation of 
an optically thin plasma, we also assumed the presence of 
photoabsorption in the source's spectrum with a column 
N H L = 5 x 10 22 cm~ 2 . Since PCA/RXTE observations 
are available after 0.7 day from the outburst onset, we 
used the temperature measured precisely by PCA in this 
outburst period. 

We clearly see from the ASM and PCA data that the 
radiation temperature during the first two days remained 
approximately constant and equal to ~ 10 keV (Fig.l). 

In the case of a strong shock wave, the matter tem- 
perature at the front is essentially determined only by its 
velocity: 

2 7 - 1 
fcTi = 2[im p D - — ——r, 
(7 + 1) 

where 7 is the adiabatic index of the interstellar gas 
and D is the shock velocity. Consequently, the shock ve- 



locity can be estimated from the measurements of the ra- 
diation temperature at the shock front: 



D = ( 7 +l) i 



fcTi 



2(7 - l)fim p 



= 28M(T 1A0keV ) 1/2 km/ Sl 



where T ljl0ke v = T hkeV /10. 

However, in our X-ray observations, we measure not 
the plasma temperature at the shock front but the flux- 
averaged temperature of the entire post-shock matter. The 
relationship between the mean temperature and the tem- 
perature at the shock front depends on the distribution of 
the pre-shock stellar wind density and the radiative cool- 
ing of the post-shock heated matter. Let us consider how 
these factors affect the relationship between the tempera- 
tures. 

7.1. Shock Motion Through a Region with a Uniform 
Stellar Wind Density 

When the piston moves with a constant velocity through a 
stellar wind with a constant density, the post-shock tem- 
perature is constant and depends weakly on time if the ra- 
diative cooling of the matter is negligible. Figure 11 (left) 
shows the post-shock temperature profiles 1.5 and 3 days 
after the explosior0 for an initial density of the medium 
n = 10 9 cm~ 3 . The radiative cooling of the matter at such 
density is negligible, but a reduction in the temperature 
near the inner boundary due to radiative losses is already 
noticeable on day 3. The time dependences of the mean 
post-shock temperature (solid line) and the temperature 
at the shock (dashed line) in this calculation are shown in 
Fig. 11 (right). The theoretical temperature at the shock 
is 10.5 keV. We see from Fig. 11 that the mean temper- 
ature within the first several days does not differ greatly 
from the temperature at the shock and, to a first approxi- 
mation, the observed temperature can be set equal to the 
temperature at the shock to estimate the shock velocity. 

For the observed temperature kT ~ 10 keV in the case 
of negligible radiative losses, we obtain the velocity of the 
forward shock in the initial expansion period (until ~2 
days from the explosion onset), D ~ 2900 km/s. 

In turn, the velocity of the strong shock ahead of 
a spherical piston in a homogeneous medium is deter- 
mined by the piston velocity as D ~ 1.1U (Sedov 1945). 
Consequently, the envelope velocity in our case is U ~ 
2630 km/s. 

7.2. Shock Motion Through a Region with a Stellar 
Wind Density Profile ~ r~ 2 

When the shock moves with a constant velocity through 
a medium with a stellar wind density profile (no = 5 x 



1 The temperature rise near the piston envelope is a numer- 
ical effect (see, e.g., Noh 1978); this is the so-called "entropy 
wake" (Samarsky and Popov 1992). 
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Fig. 10. Left top: time dependence of the mean matter temperature (solid line) and the post-shock temperature 
(dashed line) when the piston moves with a constant velocity in a medium with a stellar wind density profile n ~ sinr; 
right top: time dependence of the luminosity; left bottom: deviation of the mean matter temperature from the mean 
value; right bottom: deviation of the plasma luminosity in this calculation from the law L = At 3 expected if the 
pre-shock matter were distributed uniformly. 




Distance, 10 13 cm Time, days 

Fig. 11. Left: post-shock temperature profiles when the envelope moves with a constant velocity through a medium 
with a constant density n = 10 9 cm -3 1.5 and 3 days after the onset of envelope expansion; right: time dependences 
of the mean temperature (solid line) and the temperature at the shock (dashed line). 



10 10 (r/5.1 x 10 n )~ 2 CTO~ 3 ), the postshock temperature 
profiles 1.5 and 3 days after the onset of envelope expan- 
sion have the form shown in Fig.l2(left). The time depen- 
dences of the mean temperature and the temperature at 
the shock are shown in Fig.l2(right). 

Although the mean temperature in this case remains 
constant and we could determine the velocity of the pis- 
ton envelope using the observed mean radiation temper- 
ature, other observed quantities (e.g., the absolute lumi- 



nosity of the heated matter and its time dependence) in 
this model disagree strongly with observations. The light 
curve for this stellar wind density profile is shown in Fig. 
13; we see that the luminosity begins to decline almost im- 
mediately, while it follows from observations that it rises 
within the first ^0.75 day. 

It follows from this calculation that the stellar wind 
density profile near the white dwarf should differ from the 
law n ~ r~ 2 . 
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Distance, 10 13 cm Time, days 

Fig. 12. Left: post-shock temperature profiles when the envelope moves with a constant velocity through a medium 
with a density n ~ r~ 2 1.5 and 3 days after the explosion; right: time dependences of the mean temperature (solid 
line) and the temperature at the shock (dashed line). 
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Distance, 10 13 cm Time, days 

Fig. 14. Left: post-shock density profiles when the envelope moves with a constant velocity through a medium with 
a constant density n = 10 10 cm -3 2 and 3 days after the onset of envelope expansion. A shelf formed by the cells with 
switched-off radiative cooling is seen in the temperature profile on the third day. Right: time dependences of the mean 
temperature (solid line) and the temperature at the shock (dashed line). 




Time, days 



Fig. 13. Time dependence of the luminosity when the 
shock moves with a constant velocity through matter with 
a density profile n ~ r~ 2 . 

7.3. The Effect of Radiative Cooling on the 

Relationship Between the Mean Temperature and 
the Temperature at the Shock 

It follows from Fig. 7 that the stellar wind density near the 
white dwarf can be no ~ 10 10 cto~ 3 . At such densities, ra- 



diative cooling plays a significant role: it reduces the shock 
velocity and leads to the growth of thermal instability in 
the matter passed through the shock. 

When the envelope moves with a constant veloc- 
ity through a medium with a constant density no = 
10 10 cm -3 , the post-shock temperature profiles two and 
three days after the onset of envelope expansion have the 
form shown in Fig. 14(left). A region with switched-off 
radiative cooling in its cells (a shelf) is seen in the tem- 
perature profile on the third day near the inner bound- 
ary. Figure 14(right) shows the time dependences of the 
post-shock temperature (dashed line) and the mean tem- 
perature (solid line). It follows from this figure that in the 
case of intense radiative cooling of the post-shock matter, 
its mean temperature decreases greatly even at a constant 
piston velocity. This behavior of the mean temperature is 
in conflict with the observed one, from which it follows 
that the temperature of the emitting matter was almost 
constant for the first two days (see Fig. 1). 

Consequently, we can constrain the stellar wind den- 
sity near the white dwarf, no < 10 10 cm -3 , by assuming 
that the envelope is not accelerated. Interestingly, the de- 
rived constraint on the stellar wind density near the white 
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Fig. 15. Left: stellar wind temperature profiles when the 
shock moves with a constant velocity in a region with a 
constant density of the interstellar medium (two days after 
the explosion) and in a region where the density decreases 
as n oc r~ 2 (three days after the explosion). Right: initial 
density profile of the interstellarmedium. 




Fig. 16. Time dependences of the mean temperature 
at various stellar wind densities near the white dwarf: 
10 9 cm~ 3 (dashed line), 5 x 10 9 cto~ 3 (dotted line), and 
10 10 cto -3 (solid line). The piston velocity was kept con- 
stant (U — 2700 km/s). The boundary of the transi- 
tion from a constant density to a decreasing one is r c = 
1.9 x 10 13 cm. The dots indicated the observational data. 



dwarf agrees well with the density estimated from the 
source's luminosity (see below). 



8. Rise Phase Of The Light Curve 

In Fig. l(left), the dashed line indicates the fit to the in- 
crease m observed (3-20 keV) flux by the model oc i 3 . The 
observed rise in the source's luminosity is well described 
by the law 



(1.1 ± 0.1) x lQ 12 r s d 2 2kpcl where d 



*-2kpc 



La,-20keV 

d/2 kpc. 

It follows from our calculations that as long as the 
envelope moves with a constant velocity in a medium with 
a constant density and the radiative cooling of the heated 
matter is negligible, the post-shock temperature remains 
approximately constant. Consequently, the flux increases 
just as the hot-plasma emission measure, L oc n 2 V oc r 3 , 
which for a constant piston (and, hence, shock) velocity 
gives 

L = An 2 — (D 3 — U 3 )t 3 , 
3 



where D is the shock velocity and U is the piston velocity. 

Thus, it can be said that the observed behavior of the 
source's luminosity suggests that the shock in the binary 
system moved through matter with a constant density 
in this period. Consequently, we can estimate the mean 
density of the interstellar medium near the white dwarf. 
At temperature kT = 10 keV, the plasma emissivity is 
A3-20fceV = 3.4 x 10 -24 er<7 cm 3 s~ 1 (the value was reduced 
to the total particle number based on the adopted abun- 
dance) and the theoretical time dependence of the lumi 
nosity is L^^okeV = 93 n\ t/| 700 i 3 , where n\ is in cm 
and ?727oo = J7/2700 km/s. It thus follows that 



-3 



n {r < r c ) ~ 8.6 x 10 9 ' d 2 k P cU 27 3 ' 2 cm 3 . 

For the distance range 1.1-1.9 kpc measured by 
Barsukova ct al.(2006), the possible range of densities near 
the white dwarf is (4.8 - 8) x 10 9 cm~ 3 . 

The time elapsed from the outburst onset to the in- 
stant at which the source reaches its peak luminosity al- 
lows the distance to which the stellar wind density remains 
approximately constant to be estimated: r c = Dt pea k ~ 
1.9 x 10 13 (£ peafc /0.75da2/)A3OOO, where D 3mo = £/3000 
km/s. This value agrees with r c obtained by analyzing 
the binary's sizes (see the section above) at an orbital in- 
clination sinz ~ 0.3 — 0.4. 

9. Mean Temperature Of The Emitting Matter 
And Shock Velocity At Late Expansion Phases 

The radiation temperature begins to decrease ^2 days af- 
ter the outburst onset (Fig.l, right). The temperature can 
decrease for several reasons: 

— as a result of the shock passage into a region with 
a decreasing density of the interstellar medium at a 
constant velocity of the piston envelope; 

— as a result of the deceleration of the piston envelope. 
In this case, the shock passes to the Sedov phase (the 
regime of shock motion after an instantaneous point 
explosion). 
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Fig. 17. Time dependences of the post-shock plasma tem- 
perature (open triangles; the solid line indicates the law 
T ~ t™ 66 ) and the mean radiation temperature (open 
squares; the solid line indicates the law T ~ t~°- e ) when 
the shock moves in the Sedov phase in a medium with 
a density n oc r~ 2 . The circles indicate the observa- 
tional data. The dashed line indicates the dependence 

TkeV ~ 11 days ' 



9.1. Shock Motion Through a Region with a 
Decreasing Density 

When the envelope moves with a constant velocity and as 
the shock passes from a medium with a constant density 
to a medium with a decreasing density, the behavior of the 
emitting plasma temperature changes. Although the post- 
shock plasma temperature increases, since the shock ve- 
locity in a medium with a density profile n oc r~ 2 depends 
on the piston velocity as D = 1.1917 (Sedov 1945; Parker 
1963), the mean radiation temperature will decrease due 
to adiabatic cooling. Typical post-shock plasma temper- 
ature profiles for negligible radiative losses in this case 
are shown in Fig. 15 (left). The initial stellar wind density 
profile is shown in Fig.l5(right). 

We calculated the motion of the piston envelope with 
a constant velocity through a medium with the density 
profile shown in Fig. 15 but fixed the radius of the tran- 
sition from a constant density to a decreasing one at 
r c = 1.9 x 10 13 cm. Since the distance to the source is 
uncertain, the stellar wind density near the white dwarf 
estimated in Section "Rise Phase Of The Light Curve" 
is not the ultimate value. Therefore, below we considered 
several values of no: 10 9 , 5 x 10 9 , and 10 10 cm -3 . The time 
dependences of the temperature of the emitting matter 
derived in these calculations are shown in Fig. 16. 

We see from the figure that as a result of the piston 
motion with a constant velocity in such a stellar wind 
density profile, the mean temperature after some decrease 
again reaches a constant value at late expansion phases. 
Thus, the model with a constant envelope velocity does 
not allow the observed behavior of the temperature to 
be described at late expansion phases. Consequently, the 
piston envelope should begin to decelerate by this time, 
while the shock will pass to the Sedov phase. 



It also follows from Fig. 16 that strong radiative cool- 
ing (the curves for uq = 10 10 cm -3 ) is in conflict with 
observations even for such a density profile. 

9.2. Shock Motion in the Sedov Regime 

In the Sedov regime, the shock motion ceases to depend 
on the motion of the piston envelope: the shock velocity 
decreases with time as D ~ i~ 3 / 5 for a uniform density 
distribution of the ambient medium and as D ~ f 1 ^ 
for a decreasing (n oc r~ 2 ) density of the mcdium(Sedov 
1945). Consequently, the postshock plasma temperature 
depends on time as T ~ i~ 6 / 5 and T ~ i~ 2 / 3 , respectively, 
provided that the radiative losses of the post-shock matter 
are negligible. 

The calculated time dependences of the mean tem- 
perature of the heated matter and the temperature at 
the shock for the shock motion in the Sedov phase in a 
medium with a decreasing density, no = 5 x 10 10 (r/2.1 x 
10 11 ) _2 cm -3 , arc shown in Fig. 17. The explosion energy 
was set equal to 5.77 x 10 40 erg. The exponent in the time 
dependence of the mean temperature T ~ t~ ' 6 is close 
to both the exponent for the temperature at the shock, 
T ~ t-°- 66 , and to the observed law T ~ tr - 7 - - 6 on days 
4-10 after the explosion. Consequently, we may conclude 
that the shock in CI Cam passed to the Sedov phase on 
days 4-10 after the onset of envelope expansion. 

The explosion energy can be estimated by comparing 
the observed and model time dependences of the mean 
temperature of the emitting matter in this regime. For the 
time being, by the explosion energy we mean the kinetic 
energy of the ejected envelope, while below we will refine 
this concept. The shock velocity in the Sedov phase in a 
medium with a density hq = A\/r 2 is 

3 \ AaJ 

Here, E is the explosion energy, a = 2.1 (Book 1994), and 
the coefficient A = Aifim p ~ 10 11 n 9 r 2 13 g cm -1 , where 
rig — n /10 9 cm~ 3 and r Cjl3 = r c /10 13 cm. 

The time dependence of the temperature at the shock 

is 

T keV = 1.6 x 10- 20 [^f-) m t-£- 

It follows from the calculations presented in this sec- 
tion that the flux-averaged temperature of the post-shock 
matter is lower than that at the shock itself, but their 
ratio has the same value of 1:1.3 (here, we assume that 
the radiative losses in the plasma are negligible). The ob- 
served time dependence of the radiation temperature on 
the fourth day after the explosion can be fitted by the law 
TkeV ~ H-l^da/s 3 (indicated in Fig. 17 by the dashed line). 
Consequently, the required explosion energy is 

E ~ 2.65 x 10 42 n 9 r 2 13 erg. 
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Fig. 18. Time dependences of the mean temperature for 
the expansion of envelopes with various masses. The enve- 
lope expanded freely from the very beginning of the out- 
burst (see the text). 

10. Estimating The Mass Of The Escaped 
Envelope. Description Of The Observed 
Temperature And Luminosity 

In the previous sections, we showed the following: 

— a constant mean temperature of the emitting matter 
within the first 1-2 days after the onset of envelope 
expansion can be obtained when the piston envelope 
moves with an approximately constant velocity; 

— at late explosion phases (on days ~4-10 after the ex- 
plosion onset), the shock motion in the stellar wind 
passes to the Scdov regime. 

The shock passage to the Sedov regime takes place 
when the piston envelope transfers a substantial fraction 
of its kinetic energy to the surrounding matter and, con- 
sequently, decelerates greatly. Assuming a uniform distri- 
bution of interstellar matter with a density no , we can es- 
timate, to a first approximation, how the time ts of shock 
passage to the Sedov phase depends on the envelope mass: 

M ej ~ ^ir(DT S ) 3 fim p n 0l 

i.e. 

t s ~ 1.7MV 3 _ 6 n^ e 1 9 /3 (Z)30oo)" 1 s, 

where M e j : _6 = M e j /10 _6 Mq is the envelope mass, 
^5e9 = ^0/(5 x 10 9 )cm -3 is the unperturbed stellar wind 
density. 

However, first, it follows from our estimates obtained in 
the Section "Rise Phase of the Light Curve" that the shock 
traverses the region with a constant density within the 
first 0.75 day, the stellar matter density then drops and, 
hence, the envelope deceleration efficiency decreases; sec- 
ond, it must be known how the mean temperature behaves 
as the envelope decelerates. Therefore, we performed nu- 
merical calculations for various envelope masses in which 
its velocity was specified at the initial time; subsequently, 
it expanded freely, i.e., the envelope changed its velocity 
as a result of its interaction with the medium. The stel- 
lar wind density near the white dwarf was set equal to 



no = 5 x 10 9 c?n~ 3 and the radius of the transition from 
a constant density to a decreasing one was taken to be 
r c = 1.9 x 10 13 cm. It follows from these calculations (Fig. 
18) that an envelope with a mass M e j > 5 x 10 -6 Mo 
needed for it to expand with a constant velocity in the 
initial period of its expansion has no time to decelerate 
sufficiently at late phases, while an envelope with a smaller 
mass begins to decelerate immediately, causing the radia- 
tion temperature to drop immediately after the explosion. 
As an example, Fig. 18 shows the time dependences of the 
mean temperature for envelope masses M e j = 5 x 10~ 6 Mq 
and M ej = 1(T 6 M . 

Consequently, we may conclude that the envelope ve- 
locity in the initial period of its expansion is maintained 
not by its mass but by an additional pressure from the 
white dwarf. Most likely, one might expect this to be the 
radiation pressure of the matter that continues to burn on 
the white dwarf surface. 

For the radiation pressure to be able to effectively keep 
the envelope expansion velocity constant within the first 
1.5-2 days, the optical depth of the envelope must be much 
larger than unity. If the envelope temperature is assumed 
to be ~ 10 5 — 10 6 K, then the absorption cross section will 
be determined by the absorption of incompletely ionized 
heavy elements, which increase significantly the absorp- 
tion cross section compared to that for the scattering by 
free electrons. We cannot accurately estimate the absorp- 
tion cross section in the envelope matter, because its chem- 
ical composition and ionization fraction are unknown, but 
we can estimate the photon absorption cross section re- 
quired for an envelope with a mass M e j = 10 -6 Mq to 
become optically thin 1.5 days after the explosion. If the 
envelope expands spherically symmetrically, filling the vol- 
ume with a uniform density n, then 

M e j ~ —R n/j,m p . 
Since r ~ naR ~ 1, we obtain: 

4-7T 2 

Af cj ~ —R fj,m p . 

OCT 

Moving with the velocity U = 2700 km/s, the envelope 
will traverse a distance R ~ 3.5 x 10 13 cm within 1.5 days 
and will become optically thin if the photoabsorption cross 
section in the matter is a ~ 2 x 10 _24 cm 2 . 

10.1. Best-Fit Models for the Observed Time 

Dependence of the Radiation Temperature 

To describe best the observed time dependence of the 
mean temperature, we must assume that the mean stel- 
lar wind density around the white dwarf remains constant 
up to the distance that the shock traverses in ~2 days, 
i.e., r c = 5.2 x 10 13 Z?3ooo cm. Figure 19 shows the best- 
fit models for r c = 5.2 x 10 13 cm and two matter den- 
sities near the white dwarf: no = 10 9 cm -3 (dashed line) 
and no = 5 x 10 9 c7n~ 3 (solid line). However, without de- 
scribing these models in detail, we will immediately say 
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Fig. 19. Left: time dependences of the mean temperature when the envelope expands in a medium with r c = 5.2 x 10 13 
cm. The solid and dashed lines indicate the dependence for tlq = 5 x 10 9 cm -3 and no = 10 9 cto -3 , respectively. Right: 
time dependence of the luminosity. 



that such a value of r c is possible at low orbital inclina- 
tions (sini < 0.1). In addition, the main shortcoming of 
these models is that such a density distribution causes 
the X-ray luminosity to increase as L oc t 3 until approxi- 
mately the same time (2 days since the onset of envelope 
expansion). This is in conflict with the observational data 
(Fig. 19, right) suggesting that the luminosity reaches its 
peak ~0.75 day after the outburst onset. Therefore, be- 
low we assume that the region with an approximately con- 
stant stellar wind density around the white dwarf extends 
to r c = 1.9 x 10 13 cm; subsequently, the density falls as 
n oc r~ 2 . 

As a result of the uncertainty in the stellar wind den- 
sity no near the white dwarf and the uncertainty in the 
time of envelope motion under the action of an exter- 
nal force, several models with different envelope masses 
describe the time dependence of the mean temperature 
equally well (Fig. 20). The model parameters are given in 
the table; as the initial conditions, we specified the stellar 
wind density near the white dwarf no, the envelope mass 
M e j, and the envelope pushing time At in our calcula- 
tions. Also given in the table are the kinetic energy of the 
envelope E e j and the total energy of the post-shock stel- 
lar wind E sw at the time when the external force ceases 
to act and a theoretical estimate of the explosion energy 
via the observed temperature at the envelope deceleration 
phase E se dov In these models, by the explosion energy 
we mean the kinetic energy of the envelope and the to- 
tal energy of the post-shock matter at the time when the 
envelope ceases to be affected by the external force. The 
envelope velocity in the initial period of its expansion was 
U = 2700 km/s. 

We see from the table that the kinetic energy of the 
envelope in model 1 is lower than the total energy of the 
post-shock matter; accordingly, it is not the main source 
of the shock energy at late phases. Since a decrease in 
the envelope mass in this model will cause the mean tem- 
perature in the Sedov regime to change only slightly, the 
derived envelope mass is most likely an upper limit in this 
case. 
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Fig. 20. Time dependences of the mean temperature for 
models with various initial conditions: model 1 (solid line), 
model 2 (dotted line), and model 3 (dashed line). The 
models are described in the text. The dots indicate the 
observational data. 

It should be noted that although the shock passed to 
the deceleration phase at late expansion stages in models 
1 and 3, this is not yet the Sedov regime, when the post- 
shock temperature depends on time as ~ i~ 2 / 3 . 

We sec from Fig. 20 that there are no qualitative dif- 
ferences in the description of the time dependences of 
the mean temperature by these models. Nevertheless, a 
small envelope mass makes it difficult for the envelope 
to move under the action of an external force within 1.5 
days, because it rapidly becomes optically thin. Therefore, 
we will consider the range of possible envelope masses 
M e j ~ 10~ 7 - 10- 6 M S . 

10.2. Accelerated Motion of the Envelope 

The shock passage into the region with a decreasing den- 
sity ~0.75 day after the outburst causes a decrease in the 
mean temperature and a discrepancy between the calcu- 
lated and observed temperatures on the second day of en- 
velope expansion. In order to try to remove this discrep- 
ancy, wc performed calculations with accelerated envelope 
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Fig. 21. Time dependences of the mean temperature for 
accelerated motion of the envelope (see the text).T he dots 
indicate the observational data. 



motion. It follows from these calculations that for the tem- 
perature to be kept constant within the first two days, the 
envelope velocity must be increased by a factor of ^2 by 
the second day of its expansion, but as a result of this 
measure, the calculated mean temperature turns out to 
be much higher than the observed one at late expansion 
stages. In Fig. 21, the solid line indicates the time depen- 
dence of the mean temperature for the model with the fol- 
lowing parameters: the stellar wind density near the white 
dwarf was no = 5 x 10 9 cm _3 (Y c = 1.9 x 10 13 cm); the en- 
velope velocity was kept constant and equal to U = 2700 
km/s within the first 0.8 day; subsequently, the velocity 
increased in such a way that it was twice that at the onset 
of acceleration by the second day; after the second day, the 
envelope expanded freely, with the envelope mass after the 
onset of its free expansion having no effect on the shock 
motion. It follows from these calculations that the absence 
of envelope acceleration in our calculations is not respon- 
sible for the discrepancy in temperatures on the second 
day of expansion. 

10.3. Description of the Binary Luminosity During its 
Outburst 

The time dependences of the luminosity for models 1, 2, 
and 3 are shown in Fig. 22. The discrepancy between the 
calculated and observed luminosity at the rise phase is 
unimportant, since there is an uncertainty in the esti- 
mated stellar wind density near the white dwarf due to 
the uncertainty in the distance to the system. At the same 
time, the discrepancy in the time dependences of the lu- 
minosity at the decline phase can be explained by the fact 
that the behavior of the luminosity during the outburst is 
subjected to density variations much more strongly than 
the behavior of the radiation temperature. Therefore, in- 
complete agreement between the model and observed time 
dependences of the luminosity during the 1998 outburst 
of CI Cam is most likely related to peculiarities of the 
actual stellar wind density distribution in the binary sys- 
tem, which we do not know at present and which can- 




Time, days 



Fig. 22. Time dependences of the luminosity for various 
initial conditions: model 1 (solid line), model 2 (dotted 
line), and model 3 (dashed line). The dots indicate the ob- 
servational data. 




Time, days 

Fig. 23. Time dependences of the system's luminosity for 
a spherically symmetric stellar wind (solid line) and a stel- 
lar wind with a dense disk of constant thickness (dashed 
line). The dots indicate the observed luminosity 

not be taken into account in spherically symmetric calcu- 
lations. For example, one might expect the stellar wind 
density distribution around the B[e] star to be nonspheri- 
cally symmetric and more closely resemble a disk. In this 
case, the heated matter of the stellar wind disk will make 
the main contribution to the observed X-ray emission. A 
disk with a constant opening angle will only reduce the 
luminosity normalization but will not change the pattern 
of luminosity behavior qualitatively However, if the disk 
has a constant thickness (see, e.g., Ignace et al. 2004), 
then the decline in luminosity during the outburst will be 
faster. 

Figure 23 shows the time dependences of the system's 
luminosity for the shock motion in a medium with a spher- 
ically symmetric density distribution around the white 
dwarf and under the assumption that dense matter in the 
form of a disk with a constant thickness makes the main 
contribution to the emission at the decline phase of the 
light curve. The luminosity at the decline phase for a disk 
wind in this figure was estimated only qualitatively (the 
normalization of the light curve at the decline phase was 
chosen in such a way that the luminosity at the peak of the 
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Fig. 24. Temperature profiles downstream of the forward 
and reverse shocks. The forward shock is at r ~ 4.5 x 10 13 
cm, the reverse shock is at r ~ 3.8 x 10 13 cm, and the 
contact boundary is at r ~ 3.9 x 10 13 cm. 

light curve coincided with its value from the light curve at 
the rise phase) by assuming that the time dependence of 
the temperature did not change compared to the spheri- 
cally symmetric case. 

11. Contribution Of The Matter Heated By The 
Reverse Shock To The Observed Temperature 
And Luminosity 

The time dependence of the temperature downstream of 
the reverse shock (RS) is affected by many factors that 
are not known exactly to date, such as the density profile 
in the envelope; whether the envelope expands, having a 
finite thickness, or continuously fills the entire space up to 
the white dwarf. In this paper, we will give only qualitative 
estimates based on the decay of an arbitrary discontinu- 
ity (Riemann problem) in the plane case (Rozhdestvensky 
and Yancnko 1978). These will show that the contribution 
from the post-RS matter to the observed emission may be 
neglected, to a first approximation. Assuming the expan- 
sion of the envelope to be uniform, its density depends on 
the radius as 



M, 



>j 



'<•] 



4/37rr 3 [im p 



= 3.7 x 10 15 Af_ 6 , ej r- 3 1011 



Under certain conditions, a reverse rarefaction wave 
can be generated in the envelope at the initial time, which 
will subsequently turn into a shock as the envelope ex- 
pands. Let us estimate the distance at which the RS gen- 
eration condition is satisfied for the derived envelope ve- 
locity of 2700 km/s. The dependence of the velocity of the 
envelope at which the RS will travel into it on the ratio of 
the pressures in the envelope p e j and the ambient medium 
Po is 



U shock = (1 - h)co 



'Pej/po 



1 + h 



Pej/P0 + h 



where h = For the adopted stellar wind tempera- 

ture kTo = 1 eV, the speed of sound is cq = 16 km/s. 



For our estimates, we will take the stellar wind den- 
sity riQ = 10 10 c?7i~ 3 . For adiabatic envelope expansion, 
the pressure ratio depends on the radius as p e j/po = 
"ej,o/no(5 x 10 n /r) 37 = 3.7 x 10 5 r^ 1Qll , where n ej . Q cor- 
responds to the envelope density at r — 5 x 10 11 cm; the 
envelope and stellar wind temperatures at this radius were 
also assumed to be equal. 

Assuming that p e j/po » 1> w ^ find 



Pej/P0 



U. 



shock 



co(l-h) 



(l + ft)~6.3x 10 4 



The distance that the envelope shoul reach in or- 
der to have such a density (for the M e j = 10" 6 M Q ) is 
r = 7 x 10 11 cm , i.e., the RS is generated ~0.7 h af- 
ter the onset of expansion. As the envelope expands, the 
pressure decreases, causing the RS to strengthen. By the 
shock strength we will mean the ratio of the shock velocity 
to the speed of sound in the medium. Let us estimate the 
time when the RS will be able to heat the matter up to 1 
keV. The RS strength M r depends on the strength of the 
forward shock (FS) Mf as 



Mi 



where a = po/pej and (3 



aMj + [3, 



h Pe 



-i>n 



if the adiabatic 



l+h Pe 

indices of the matter in the ambient medium and the 
envelope are equal; M r = Mf = ^ 



it = =- ~ 187.5, 
, e] — \D\ — U\\ is the RS velocity in the envelope matter 
frame of reference. 

For the RS to be able to heat the matter to 1 keV, 



its velocity relative to the envelope must be D e 



915 



km/s. It follows from the relationship between the shock 
strengths that 



Po 

Pej 



The pressure in the envelope can be expressed in terms 
of p e j = 7 _1 jj,m p n e jc1j. Thus , for the RS to be able to 
heat the matter to 1 keV, the density in the envelope must 
be n e j = 1.1 x 10 n cm~ 3 . Consequently, the envelope must 
expand to the distance r = 1.6 x 10 13 cm, i.e., ^0.7 day 
should be elapsed since the onset of expansion. Let us es- 
timate the transparency of the post- RS matter on 0.7 
day after the onset of expansion if the main mechanism 
responsible for the matter opacity is Thomson scattering: 
r - Axn ej a T Ar - 4x 1.1 x 10 11 x 6.65 x 10~ 25 x 3 x 10 12 ~ 
0.9. For the absorption cross section a = 2 x 10 _24 cm 2 
obtained above, r ~ 2.6. As long as the heated post- 
RS region is optically thick, the entire envelope emits the 
Eddington luminosity and the temperature of the emer- 
gent radiation is very low, ~ 3 eV. 

The envelope becomes optically thin as it expands. 
To demonstrate the behavior of the post-RS temperature, 
Fig. 24 shows an example of the temperature profile in 
the case where the envelope becomes optically thin almost 
immediately and cools down through the radiation of an 
optically thin plasma (the FS is at r ~ 4.5 x 10 13 cm, the 
RS is at r ~ 3.8 x 10 13 cm, and the contact boundary is 
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Table 1. Parameters of the best- fit models for the observational data 



Model 
number 


At, days 


n , 10 9 cm -3 


M ej , M Q 


£ ej ,10 43 erg 


E sw ,10 43 erg 


E sedov ,10 i3 erg 


1 


1.5 


5 


2.85 x 1(T' 


2 


3.3 


4.8 


2 


1 


5 


1(T 6 


8 


1.7 


4.8 


3 


1.5 


8 


i(r 6 


8 


5 


7.7 



at r ~ 3.9 x 10 13 cm). It follows from this profile that as 
long as the envelope density was high (the region formed 
near the contact discontinuity), the temperature of the 
post-RS matter was low, since it effectively cooled down 
through radiation, but the post-RS matter temperature 
began to increase as the envelope expanded and became 
rarefied. We calculated the mean temperature of the radia- 
tion coming only from the FS-heated matter and the mean 
temperature of the entire emitting matter for this temper- 
ature profile. As we suggested, the contribution from the 
post-RS matter does not affect the mean temperature and, 
hence, it will not affect the luminosity in the observed 3-20 
keV energy band. 

Thus, our preliminary estimates show that, to a first 
approximation, the contribution of the radiation from the 
RS-heated matter to the observed 3-20 keV luminosity 
and the mean temperature is small compared to that of 
the radiation from the FS-heatcd matter. 



12. Conclusions 

Using analytical estimations and numerical simulations, 
we were able to describe the intense X-ray outburst in 
the binary system CI Cam observed in 1998 in terms of 
the model of the interaction between the expanding enve- 
lope ejected as the result of a thermonuclear explosion on 
the white dwarf surface (classical nova explosion) and the 
stellar wind from the optical component. In particular, we 
showed the following: 

— The mean radiation temperature of the stellar wind 
matter passed through the shock differs significantly 
from that of the matter immediately behind the shock 
front. Therefore, numerical calculations are needed to 
properly determine the classical nova envelope expan- 
sion parameters. 

— According to our preliminary estimates, the contribu- 
tion from the reverse shock wave heated matter to the 
observed mean temperature and the 3-20 keV luminos- 
ity is smaller than that from the forword shock wave 
heated matter. 

— The radiation temperature measured at early envelope 
expansion phases of the classical nova allowed its ex- 
pansion velocity to be estimated: ~2700 km/s. 

— In our model, the envelope is ejected from the white 
dwarf as the result of explosive thermonuclear burning. 
It has an expansion velocity of ^2700 km/s already on 
0.1-0.5 day after the explosion onset and, being opti- 
cally thick, moves under the radiation pressure from 



the white dwarf with a constant velocity for the first 
~1-1.5 days. Thus, we have been able to measure the 
envelope velocity almost immediately after the explo- 
sion for the first time. 

— Subsequently, the envelope probably becomes optically 
thin and decelerates while interacting with the stellar 
wind from the optical component. The shock in the 
stellar wind passes to the Scdov phase with the time 
dependence of the mean temperature T ~ £-0.7-o.6_ 

— By comparing the observed rise in luminosity with 
the theoretical dependence, we estimated the stellar 
wind density near the white dwarf, no(r < r c ) ~ 

8.6 x 10 9 d2fc pc L^27ou 2cm ~ 3 - I n the simplest model of 
the stellar wind density distribution, this corresponds 
to a stellar mass loss rate in the stellar wind of ~ 
(1 — 2) x 10~ 6 AIq /yr. This estimate agrees well with 
the values obtained by Robinson ct al. (2002). 

— The observed time dependence of the temperature of 
the emitting matter at late envelope expansion stages 
allowed the mass of the ejected envelope to be con- 
strained in our model, ~ 10 — 10 _6 M Q . 
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